Application of the density dependent hadron field theory to 

neutron star matter 



F. Hofmann, C. M. Keil, H. Lenske 
Institut fur Theoretische Physik, Universitat Giefien, Heinrich- Buff- Ring 16, 35392 Giefien, 

Germany 
(February 8, 2008) 



Abstract 



The density dependent hadron field (DDRH) theory, previously applied to 
isospin nuclei and hypernuclei is used to describe /3-stable matter and neutron 
stars under consideration of the complete baryon octet. The meson-hyperon 
vertices are derived from Dirac-Brueckner calculations of nuclear matter and 
extended to hyperons. We examine properties of density dependent inter- 
actions derived from the Bonn A and from the Groningen NN potential as 
well as phenomenological interactions. The consistent treatment of the den- 
sity dependence introduces rearrangement terms in the expression for the 
baryon chemical potential. This leads to a more complex condition for the j3- 
equilibrium compared to standard relativistic mean field (RMF) approaches. 
We find a strong dependence of the equation of state and the particle distri- 
bution on the choice of the vertex density dependence. Results for neutron 
star masses and radii are presented. We find a good agreement with other 
models for the maximum mass. Radii are smaller compared to RMF mod- 
els and indicate a closer agreement with results of non-relativistic Brueckner 
calculations. 



PACS number(s): 26.60.+C, 21.65.+f, 21.30.Fe, 97.60.Jd 
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I. INTRODUCTION 



The nuclear equation of state (EoS) is the fundamental input for the calculation of 
neutron star properties. The internal structure of a neutron star ranges from sub-nuclear 
densities at the surface to a few times the normal nuclear matter density po i n its core. 
Therefore, a detailed knowlegde of the EoS over a wide range of densities is required. Of 
particular importance is the behavior of the EoS for densities p ^> po since it primarily 
determines the maximum mass of the star. Comparision with experimentally observed data 
for neutron star masses and radii might allow us to examine the reliability of different models 
for the EoS at high densities. 

The theoretical determination of such an EoS is a very hard task since there is only 
accurate knowledge around p = po where matter in the weak equilibrium consists mainly 
of nucleons and leptons. At higher densities of 2 — 3 po one expects strange baryons to 
appear as new hadronic degrees of freedom. The extrapolation of the EoS to high densities 
and its implication on neutron star properties has been examined in a variety of models, 
most notably in the framework of phenomenological non-relativistic potential models [|TJ, 
relativistic mean-field approaches (RMF) 0-[5[] and non-relativistic Brueckner-Hartree-Fock 
(BHF) calculations ||[7|]. In addition there have been calculations employing chiral effective 
Lagrangians || or the quark meson coupling model |J. From a microscopic point of view 
Brueckner calculations using realistic nucleon-nucleon (NN) and hyperon-nucleon (YN) po- 
tentials as input have to be preferred. Since the density inside neutrons stars is extremely 
high - the Fermi momenta and the baryon effective mass are in the order of 500 MeV - 
the use of non-relativistic models might be problematic and one should prefer a relativistic 
description. Due to technical difficulties, relativistic (Dirac) BHF calculations are restricted 
to asymmetric nuclear matter and the inclusion of hyperons at higher densities is at present 
not feasible. Therefore, only DB calculations for neutron star matter with nucleons have 
been performed [10jll.|. 



RMF theory fl"2| , [T3"|j allows to incorporate easily and consistently an enlarged set of 



hadronic degrees of freedom and has been succesfully applied to hypernuclei ||14|-|T7H and 
neutron stars @|§,f|]. Being a phenomenological model with its parameters usually adjusted 
to the properties of finite nuclei and their limited density range around p , an extrapolation 
to higher densities has some uncertainties. Parameterizations with an excellent describtion 
of finite nuclei turn out to be unstable at higher densities being mainly caused by divergent 



scalar self-interaction terms [18]. Nucleons aquire negative effective masses and the equation 



of state turns out to be much stiffer compared to Dirac-Brueckner calculations ||. One 



solution is the introduction of a quartic vector self-interaction term |fL9| , pO|l . Also, calculations 
in relativistic Hartree-Fock approximation with coupling constants fitted to the EoS of DB 
calculations were succesfully applied to neutron star matter 0,15 . 

In this paper we employ results from DB calculations by parameterizing the DB self- 
energies in terms of density dependent coupling functionals and apply them in Hartree 
approximation to strange and neutron star matter. This approach, the density dependent 
hadron field (DDRH) field theory, has been introduced in [^3|,f24[] and applied to stable 
[25-p7| and exotic nuclei f28fl . Taking into account information from realistic NN potentials 
and a much wider density range (0.5 — 3 po) we expect our extrapolation to high densities 
to be more stable than purely phenomenological RMF models. 
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An important point is the treatment of the hyperons. In the DDRH model they can be 
treated in almost the same manner as in standard RMF theory, allowing an easy extension 
of the model to the SU(3)/ octet baryons as was shown in [E{|. Since we assume in general 
density dependent hyperon-meson vertices we should in principle derive their parameteriza- 
tion from DB calculations that include hyperons. Since such calculations are not available 
we extrapolate the density dependence and the strength of the vertices from DB calculations 
of nuclear matter and from hypernuclear data. As discussed in [T2^] this provides at least 
qualitative information about DB results of strange matter. This approach is discussed in 
Sec. [TI| where also a short review of the theoretical model is given. In Sec. [TTT] we apply the 
derived density dependent hyperon-nucleon interaction to iVA matter and to neutron star 
matter in /3-equilibrium. We discuss the influence of the density dependence on the chemical 
potential and present results for the equation of state for various density dependent interac- 
tions. Neutron star matter compositions for different models are presented. In Sec. [IV| the 
influence of the calculated EoS on mass-radius relations of neutron stars is investigated and 
compared with previous results of other groups. The paper closes in Sec. |V| with a summary 
and conclusions. 



II. DENSITY DEPENDENT HADRON FIELD THEORY WITH HYPERONS 

A. The Model Lagrangian 

The model Lagrangian closely follows in structure the one of relativistic mean- field 
(RMF) theory [12|,[L|. The important difference in the density dependent relativistic hadron 
field theory (DDRH) corresponds to replacing the constant meson-baryon vertices of the 
RMF model by functionals T(p) of Lorenz-scalar bilinear forms p(^F,^p) of the baryon 
field operators. This step is necessary to retain thermodynamical consistency and energy- 
momentum conservation for this extended model and has been discussed in detail in 



In the meson-exchange particle sector we have also included the scalar isovector meson S in 
contrast to previous formulations since it is expected to be important in very asymmetric 
systems. 

In this work the extension to hypernuclei as introduced in [29] is used. The l/2 + baryon 



octet is taken into account including the S = —1 (A, E) and S = — 2 (H) hyperons. Besides 
the standard set of non- strange mesons we also include the hidden-strangeness meson fields 
a s (scalar, m as = 975 MeV) and (vector, m = 1020 MeV) [0. 
Defining the flavour spinor typ 

^ F = (^7V,^A,^E,^ S ) T (1) 

which is composed of the isospin multiplets 
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the Lagrangian is written as 



£ — £b + £*m + £-int 
C B = ^ F [i 7/t ^ - M F ] V F (3) 



2 ■ < 



\ E (~^ ( ^ - m\AfA*A (4) 

K=UJ,p,(f> 

C int = ^pT^p^ - VfT^fA^ + 

^ F f - VfT^fA^. (5) 

and £m are the free baryonic and mesonic Lagrangians, respectively. Baryon-meson 
interactions are described by C int that includes the vertex functionals T(\l/ F , \I/ F ). The diag- 
onal matrix M contains the free-space baryon masses. The Lagrangian and the interaction 
have to be symmetric in flavor space, since SU (3) /-flavor exchanging mesons are not con- 
sidered. This is obtained by defining the vertices as 

(r a ) BBI =T aB 5 BB ,. (6) 

The indices a = a, u, 5, p, a s , and B = N, A, E, S denote all mesons and baryon multiplets. 
The most general ansatz for the DDRH vertices, allowing to treat the density dependence 
of each vertex independently, is 



aB 



(p qB (*f,Vf)) (7) 



where p aB is a Lorentz-scalar combination of the baryon field operators. As discussed in 
the introduction, the strength and the intrinsic density dependence of the vertices have 
to be deduced from microscopic calculations. The mapping of the DBHF self-energies to 
infinite nuclear matter DDRH vertices is done in the local density approximation (LDA), 



e.g. |23lj31| , |32|1 , and has been thoroughly discussed in |28|j29| . In |?8| a momentum correction 
of the self-energies was introduced to improve the applicability of DB calculations to the 
relativistic Hartree approximation. DB self-energies derived in asymmetric nuclear matter 



32[ were parameterized and special care was taken to reproduce the DB binding energy 
over the complete asymmetry range from isospin symmetric nuclear matter to pure neutron 
matter. 

The standard choice is to let the density operator p depend on the baryon vector current 
= ^7^^/. This so called vector density dependence (VDD) leads to quite satisfactory 
results for finite nuclei [2lJ and hypernuclei and is a natural choice for the parame- 



terization of the DB self-energies as was discussed in [|28| . A straightforward extension to 
meson-hyperon vertices is the ansatz [E9| 



p aB [^ F , ^f] = ^fKbIh^f- (8) 

Different choices of the matrix B^ B and their physical significance will be discussed in the 
next section. Taking the variational derivative of the Lagrangian 
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5C 



int 



dCi n t dCi n t Sp a B 



(9) 



leads to the usual meson field equations known from RMF theory but replacing the constant 
couplings g a B by density dependent vertices r aB (p). In the baryon field equations additional 
rearrangement contributions appear that have their origin in the second term of Eq. (Pf), 



7^ (id" - S M(0) - t^ r) ) - (M - E' 



(0) 



0. 



(10) 



The scalar self-energy has the standard form, while the vector self-energy 

= f>(o) + £n(r) inc i uc ies rearrangement contributions £^( r ) introduced by the medium- 
dependence of the vertices. 



E^ 0) = f W A^ + f TA ip) v + f 6 A w " 



|]A«(r) _ 



dp a 



aB 



(11) 

(12) 
(13) 



The importance of the rearrangement energies £^ r ) and their physical origin has been 



discussed in detail in |25|,|33|]. Their explicit form and their impact on neutron star matter 
will be discussed in the next sections. 



B. Choice of the density dependence 

The density dependence of the hyperon-meson vertex functionals should in principle be 
derived directly from DB self-energies of strange matter performed with the complete baryon 
octet. However, such a full scale calculation is not available and hardly feasible under present 
conditions. Non-relativistic Brueckner calculations with microscopic interactions have been 



performed for strange matter by several groups p^- pq] but, obviously, their results can not 
be used in relativistic calculations because of the different structure of the non-relativistic 
single particle potentials. Therefore, we choose a semiempirical approach by relating the 
microscopic density dependence of the meson-nucleon vertices derived from DB calculations 
of pure isospin matter (only p, n) to the density dependence of the hyperon- nucleon and 
hyperon- hyperon interaction. 

In [^9| we have shown by inspecting the properties of the Dirac- Brueckner interaction 
in strange matter that density dependent nucleon and hyperon dynamics can be related to 
each other to a good approximation by scaling laws. The main outcome is that the hyperon 
and nucleon self-energies and vertices are related to each other by the ratio of the free space 
coupling constants g a s- in leading order Hartree approximation one finds the relation 

R aY = ^ - ^ - K (14) 

Since this relation is strictly valid only for symmetric hypermatter with the same content 
of nucleons and hyperons p9f it is not obvious how the density dependence of the hyperon 
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vertices has to be related to the full baryonic density for other matter compositions. An 
obvious and rather simple approach is to let all vertices depend on the total baryon density, 
denoted by px- From now on we will denote this choice as model 1. For the hyperons this 
corresponds to defining the coupling functionals as 

TaS = RaT^ aN{pT)i 

f a Z = -Ra=f Q jv(p T ). (15) 

The parameterization of the density dependence is taken from the nucleon- vertices r Q 7v(p7v) 
derived from DB calculations of nuclear matter where p^ = p?- By relating the density px 
to the total baryon vector current 



Pt = yjjp^- (16) 

a Lorentz-invariant expression is obtained which is equivalent to choosing the matrix B^ B 
in Eq. (||) as 

Bf = B^ = V* ia g(M,M), (17) 

where w M is a four velocity with u^u^ = 1. This choice leads to rearrangement self-energies 
of the form 



B t \ ORB' OPB' 



Op B > Op B ' 

- ^*.JH,) V (18) 

OpB' OpB' J 

where the scalar densities are defined as p s B = ^Ib^b- Since the density dependence of all 
vertices was chosen to depend on the same argument pr the summation has to be performed 
over all B' leading to identical rearrangement self-energies for all baryon multiplets. Also, 
in a system consisting of different fractions of hyperons and nucleons, the scaling relation 
for the vertices is fulfilled for all densities and independent of the strangeness fraction. 
In Ref. |29| A-hypernuclei were considered and the medium modification of the vertices 



was chosen to depend only on the density of the surrounding baryons of same particle type, 
e.g. T a A{PA)- The parameterization of the density dependence was taken from the nucleon- 
vertices T a N{PN) derived from DB calculations of nuclear matter. We will denote this choice 
as model 2. Extending this model to the complete baryon octet, we let the vertices only 
depend on the density of baryons of the same multiplet 

FqA = RaA^aN(P\), 
TqS = Ra^aN^Pi:), 

r aH = RastaN(ps), (19) 
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where the densities are the baryon vector currents of the corresponding baryon species 



In this model the matrix B^ B in Eq. (|8|) is defined as 



Pb = Jtiji*. (20) 



Bf ee B B = ^diag (S NB , 5 AB , 6* B , S SB ) (21) 
and the rearrangement self-energies now differ between the baryon multiplets 




BY 



~A^j B -^-^ s p B )u^. (22) 



dp 



B 




This separation of the density dependence also corresponds in first order to the outcome of 
DB considerations and leads to satisfying results for A-hypernuclei. It takes into account 
that in leading order the medium dependence of the vertices is only caused by Pauli-blocked 
intermediate states of baryons of the same multiplet. On the other hand, this separation 
leads to a variation of the relative strength of the vertices in strongly asymmetric systems, 
e.g. neutron stars, since the vertices depend on different arguments p B . This behavior 
is different compared to standard RMF calculations where the ratio of the vertices stays 
constant over the whole density range and is independent of the strangeness content. 

The two vertex models constitute limiting cases where model 1 (dependence on px) is 
probably most realistic for systems at high total baryon densities. Modell 2 (dependence 
on pb), on the other hand, can be expected to be most realistic for low hyperon densities 
as found in single A-hypernuclei. We assume that a realistic density dependence should be 
composed as a mixture of the total density and the baryon multiplet density 

f a s = t aB (pB,PT) (23) 

Unfortunately, DB calculations from which such a functional dependence could be ex- 
tracted are not yet available. Our choice, taking into account both extremes, allows us 
to examine thoroughly the properties of such a parameterization. Furthermore, for pure 
systems, consisting only of e.g. A, both models are identical. 



C. The vertex scaling factors 



In RMF theory the phenomenological hyperon and nucleon vertices can be related to 
each other by simple scaling factors R a y, e.g. g a \ = R Q Ag a N- As shown in Section [II B 



this relation also holds in the DDRH model. A widely used approach is to determine these 
vertex scaling factors from SU(6) symmetry relations of the quark model [37]. For the vector 
mesons ideal mixing is assumed and one finds BO 
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2T 



2T 



4>A 



2T 



0S 



<j>N 



2r p 7v, 

2i> 
o. 



r pA = o 

2\/2 



r 



uiN 



(24) 



(25) 



The scalar coupling constants are fitted to hypernuclear properties [p^-|16| and chosen to 
give hyperon potentials in saturated nuclear matter that are compatible with experimental 
results for the single particle spectra of hypernuclei. 



U 



(AO 



U- 



(AT) 



-30 MeV, U-. 



(AT) 



-28MeV. 



(26) 



The values of the A, E and H potentials were chosen in accordance with Ref. PJSOf . Based 
on the analysis of S~ atomic data the real part of the optical potential was found to be 
negative |38| . A recent analysis, however, indicates that the isoscalar potential changes sign 
in the nuclear interior and becomes repulsive [3J|. Therefore, we also performed calculations 
with a potential depth of = +30 MeV and found that the £ does not appear in the 
composition of neutron star matter. For this reason we restrict our discussion to the negative 
potential value. The H nuclear interaction also exhibits large uncertainties with potential 
depths ranging from —14 MeV to —28 MeV in some investigations M| . 

In DDRH theory the non-relativistic single particle potentials of hyperons in nuclear 
matter are obtained as 



U- 



(AT) 



,0(0) 



+ E 



o(0 

Y 



,s(0) 

jy . 



(27) 



Choosing as fixed by phenomenology this relation introduces a constraint on the cou- 
plings T aY and T^y- Due to the different density dependence of model 1 and 2 one finds 
different values for the scaling factors R a y even though they were adjusted to the same 
potential depth. In model 1 the hyperon vertices have to be evaluated at saturation density 
Pt = Po even though the density of the hyperons is py -C pr- In addition, the rearrange- 
ment contributions of the nucleons add to the hyperon potential. In contrast, in model 2 
the vertices have to be calculated at vanishing hyperon density and rearrangement does not 
contribute to the hyperon potential. 

We examined the behavior of three different density dependent interactions. The Bonn 



A parameter set |3T[ parameterizes the density dependence of DB Brueckner calculations in 
symmetric nuclear matter with the Bonn A NN potential It only includes density 



dependent vertices for the a and the uj meson while the p meson coupling strength was chosen 
as a constant. Calculations for single A-hypernuclei have been performed successfully with 
this parameterization |[29|| . The Groningen parameter set [2^j was fitted to DB Brueckner 
calculations in asymmetric nuclear matter with the Groningen NN potential p2|,43| , |44 |. In 
addition to the isoscalar channel it includes a density dependence in the isovector part of the 
interaction, parameterized by the p and the 5 nucleon-meson vertices. The third parameter 
set is a phenomenological interaction whose density dependence (cr, uj and p meson) has 
been adjusted by a fit to nuclear matter properties and finite nuclei j|5| . It will be denoted 
from now on as DD. The quality of its description of finite nuclei is comparable to RMF 
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parameter sets including nonlinear meson interaction terms and its nuclear matter properties 
at higher densities are in accordance with DB calculations. Results for the sccaling factors 
are presented in the right column of Tab. |. From now on, the set of scaling factors R a Y 
that was determined from the SU(6) symmetry relations of the quark model will be denoted 
as Rq . 

It is our aim to start from DB theory and use microscopic interactions as input for our 
calculations. Since these results include highly nonlinear and nonperturbative correlation 
effects, the strict use of the quark model reduction of R^k = 2/3 might be questionable. 
Therefore, we examine in addition the properties of microscopically derived scaling factors. 
For the Bonn A potential an extension to the free A A systems exist [j4"G| , |4"?[| , but DB calcu- 
lations are pending. A free space scalar vertex scaling factor R a \ = 0.49 was extracted from 
A A T-matrix results for a sharp a meson mass m a = 550 MeV [[IT]]. Following Eq. ([14]) we 
can apply this value also to the in- medium vertices. 

Since microscopic values for R u ^ are not available from Ref. j|7J we use it as a phe- 
nomenological parameter and adjust it to the potential depth . In this analysis we are 
restricted to the A because calculations for AS and AS systems were not yet performed. 
For the Bonn A potential we find a value of R w \ = 0.569 in model 2 in close agreement 
with [29| where a value of -R^a = 0.553 was found by a fit a A single-particle energies. In a 
constant coupling RMF model a relative u coupling of R^a = 0.512 was found for the 
same value of R a \. This is again in close agreement with our corresponding values of model 
1, e.g. R^a = 0.510 for the phenomenological density dependence. Results for the vector 
meson scaling factors R^a for model 1 and 2 for the different interactions are shown in the 
left column of Tab. |. We will denote this set of semi- microscopic scaling factors by R m . 

In [fl"6"| , |4"?f the scalar meson channels were described by the correlated exchange of pion 
and kaon pairs. Therefore, the scalar coupling also includes a relevant admixtures of the a s 
field in the A coupling. Adjusting the vector coupling we also implicitly include contributions 
from the <fi field. Naive quark counting suggests that even for a pure A-A interaction the 
strange mesons contribute only about 10% of the interaction strength. For this reason we 
neglect the explicit contributions from the <j s and meson in our calculations. In addition, 
it should be noted that the contribution from these hidden- strangeness mesons can not be 
fixed by experimental data, introducing additional ambiguities in the model as seen in [40 



In the next section we will show that there are already notable differences between the 
different models without these additional fields. The notable deviations of R a A and R^a 
to the quark model value of 2/3 originate in higher order nonlinear contributions from the 
dynamically generated a and a s exchange channels P6| , |4T| , the explicit SU(3)j symmetry 
breaking and the u- <fi octet-singlet mixing on the fundamental strong interaction level. 



III. MEAN-FIELD DESCRIPTION OF STRANGE AND /3-STABLE MATTER 

A. Properties of strange matter 

We now apply the extended DDRH model of Sec. [TJ to A matter and /3-stable matter. A 
solvable scheme is obtained in mean-field theory which amounts to taking the expectation 
values with respect to the Hartree ground state of the baryon-meson vertices and the meson 
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fields. The vertices then reduce to density dependent functions r Q #(p) of the baryon density 
p. In static infinite matter all derivative terms vanish and the solution of the fields can be 
expressed analytically |T5[ . 

In weak /^-equilibrium electrons e~ and myons pT have to be included in the Lagrangian 
due to the weak decay channel, 

v|> L = (V> e ,W) T (28) 



Cl = [z7a»0" - M L ] * l . (29) 

In mean-field theory the energy density e and the pressure P are given by the ground state 
expectation values of the energy- momentum tensor T^ u . One finds 

e = (T">) = •£ \ [3E F<Pl + m*p«] (30) 



b z 



13 i 



p = tt E<n = E t - m *^l ( 31 ) 



^ t=l t=6,J ^ 

+ E^[prf 0) -P^ + 2p^ 

where the indices 6 = p, n, A, S~, S°, S + , S° and / = e~ , p~ run over all baryons and 
leptons, respectively. The self-energies are the mean-field expectation values of the quantities 
defined in Eqs. (|Tl"l) -(|l3l), Ep b = y k 2 Fb + m^ 2 is the energy of the particle, kp is the Fermi 
momentum and ml = Mb — the effective mass. 

The chemical potential in a system with fixed baryon number is defined as 

In contrast to constant coupling RMF models the rearrangement energy appears in the 
above relation which is mandatory for thermodynamical consistency. This can be easily 
verified by applying the Hugenholtz-van Hove theorem JEJ that relates energy and pressure 
of a particle to its chemical potential and noting that the rearrangement contribution also 
appears in the expression for the pressure. 

In weak /3-equilibrium the chemical potentials of all particles are related to each other 

by 

Pi = hp n - qiPe (33) 

which imposes charge and baryon number conservation. Here, the index i denotes all par- 
ticle species (baryons and leptons) and 6j and qi are the corresponding baryon number and 
electrical charge. The equilibrium problem is solved by determining the meson fields from 
the field equations of the form 
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b 



(34) 



where p depends on the choice of model 1 or 2. The densities p b are the (iso)scalar or 
(iso)vector densities and are determined by the Lorentz and isospin structure of the corre- 
sponding meson vertex. As an additional constraint the total baryon density p is fixed and 
charge neutrality is imposed for neutron star matter 

P = h iP^ Pc = Yl QiPi = 0- (35) 

i i 

This allows the determination of the electron p e and baryon p n chemical potentials in (3- 
equilibrium. In contrast to standard RMF calculations one can not simply extract the 
Fermi momentum kp. and thus the density pi = k Fi /3ir 2 of each particle for a fixed chemical 
potential from Eqs. (|32|) and (|33|) . Since the self- energies depend by virtue of the vertices 
on the densities of the baryons a coupled set of equations has to be solved self-consistently 
for the density of every baryon. In model 1 the coupling strength is known a priori from the 
given total baryon density but the rearrangement term in Eq. (^) can only be calculated 
from Eq. (|i~8D if the composition of the neutron star matter is known. For model 2 the 
calculation gets even more complex since in addition the strength of the vertices changes 
with the particle ratios. Furthermore, the rearrangement self-energies of Eq. ([22|) differ for 
each baryon multiplet and are strongly affected by the density dependence of the vertices 
and the density of each particle species. Therefore, one has to include all particles in the 
calculation and to check for their appearance such that the condition of chemical equilibrium 
is retained. This leads to a highly nonlinear set of equations where Eqs. (|32|)-(|3"5"|) have to 
be solved self-consistently and the calculation is far more involved than in standard RMF 
theory. 



B. A matter 

We first study the properties of strange matter for our different choices of the density 
dependence of the hyperon-meson vertices. This is done by calculating the equation of state 
for symmetric nuclear matter (p p = p n = \pn) with a fixed admixture fs of A-hyperons. 
The strangeness fraction for AN matter is defined as 

, Ps Pa / Qft s 

fs = — = : • (36) 

Pt Pn + Pa 



It was pointed out in (40| that AN matter will not be the lowest energetic state of strange 
matter at higher densities and for higher strangeness fractions. A full calculation ensuring 
chemical equilibrium has to include also the E and the S since the A can be converted via 
non-mesonic decay-channels into these hyperons. But for a systematic study of the effects 
of the different choices of the scaling factors R a \ and R u x for models 1 and 2 we have to 



restrict ourselves to the admixture of A-hyperons since, as discussed in Sec. [HQ , these are 
the only hyperons with microscopic derivations for the scaling factors. 

In Fig. [l| we show results for the equation of state calculated with the Bonn A parameter- 
ization for different strangeness fractions. Comparing within model 1 and 2 the differences 
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between the set of scaling factors R m from microscopic considerations and the set of quark 
model scaling factors R q , one realizes that the equation of state is slightly stiffer in the latter 
case and that the difference increases with higher densities and higher strangeness fractions. 
The reason for this is that the reduced strength of the A vertices in the microscopic case 
decreases the potential energy of the hyperons. The softer EoS is also in agreement with 
the fact that the additional hyperon-hyperon interaction, caused by the a s and <ft mesons 
and assumed to be strongly attractive at the considered densities JJfl], is at least in part 
implicitly included. At first sight this result might be surprising since in both cases 
was adjusted to the same value. But this agreement only exists at saturation density and 
for very small strangeness fractions. Since the ratio of scalar and vector couplings changes 
with the density due to their medium-dependence the equations of state differ when leav- 
ing the equilibrium density. Strictly, identical EoS were only recovered if scalar and vector 
potentials scaled proportionally to each other. This is also not the case in a RMF model 
with constant couplings since the scalar and vector potentials scale differently (partly due 
to nonlinear interaction terms) but the effect is less pronounced. 

For model 1 the minimum of the EoS is shifted to higher densities and bound stronger 
than in model 2 where the shift of the minima is practically negligible. Overall, the EoS 
is much softer in model 1 at lower densities and both models approach each other again at 
higher densities. Here, the behavior of model 2 is in closer agreement with non-relativistic 
Brueckner calculations where the position of the minimum seems to be shifted to lower 



densities [[35] or is relatively independent of the strangeness fraction [|34],[3(|. However, a 
comparison with non-relativistic results is difficult since these do not reproduce the correct 
saturation density of symmetric nuclear matter fjQfl . This behavior is closer examined in 
Fig. ^] where we display the saturation density and the binding energy at saturation density 
as a function of the strangeness fraction. The Groningen parameter set exhibits a less pro- 
nounced difference between the different models mainly due to its weaker density dependence 



at low densities |28| . The differences between model 1 and 2 are always larger than between 
the different choices R q and R m of the scaling factors. At high strangeness fractions model 
1 and 2 seem to approach each other again. This is obvious since their density dependence 
is identical in the limit of pure systems. However, the binding energy of pure A matter is 
not identical for model 1 and 2 as can be seen from Fig. ^. The reason is that the scaling 
factors have been adjusted at fs = and differ slightly due to the behavior of the vertices at 
low densities and the different choices of the density dependence. This is seen by comparing 
the scaling factors of model 1 and 2 in Tab. |. Again, the difference is most pronounced for 
the Bonn A potential. Results for the DD parameter set are not shown but closely resemble 
the Groningen results. 

This distinct difference is mainly caused by the strong density dependence of the vertices 
of the Bonn A parameterization at low and very high densities. The reason for this behavior 
is not the Bonn A potential itself but the polynomial function in fcp that was used for 
the parametrization of the self-energies. Even though this function is very accurate around 
saturation density and well suited for the calculation of finite nuclei |3l[] it leads to some 
uncertainties when extrapolated to very low densities. Since in model 2 the hyperon vertizes 
are adjusted at low densities, where due to numerical difficulties Brueckner calculations 
are problematic, the extrapolation to higher densities might involve uncertainties. The 
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extrapolation of model 1 is safe since the hyperon vertices are evaluated around saturation 
density where they are well defined. The Groningen and the DD parameter sets avoid this 
problem by choosing a rational function in p to fit the self-energies leading to a stable 
extrapolation to low and very high densities. 



C. /3-stable matter 

In the following we study the composition of the equation of state for neutron star matter 
as discussed in Sec. [Ill A| . We first examine the case where the vertices depend on the total 



density (model 1) and consider only A-hyperons. This allows us to compare the results with 



Sec. J I IB]. In Fig. £| the pressure versus the energy density is shown for the three different 
interactions. The DD parameter set has the softest EoS, i.e. the lowest pressure at a given 
e, closely followed by the Groningen parameter set while Bonn A is much stiffer at high 
energy densities. The curves calculated with the scaling factor set R q (upper lines) are 
much stiffer than the curves using the semi- microscopic set R m (lower lines). This difference 
is larger than expected from Figs. [TJ and 0. This is mainly due to the reduced repulsive 
A- vector potential for set R m . The appearance of the A sets in at slightly lower densities 
which introduces an additional softening of the EoS by a reduction of the kinetic pressure. 
Remarkably, the differences between the two choices of the scaling factors can be larger than 
the differences between two parameter sets indicating that a closer examination of the scaling 
factors R Q \ is necessary. The properties at high energy densities are mainly determined by 
the density dependence of the vertices. Here, the Groningen and DD parameter set follow 
closely Brueckner calculations. The high density extrapolation of the Bonn A parameter set 
leads to an increase of the coupling strength and the contribution from the repulsive u field 
exceeds the attraction from the a field causing the stiffer EoS. 

For comparison with constant coupling RMF calculations we also performed calcula- 
tions including the complete baryon octet. The scaling factors R a y were determined in 
the standard approach from SU(6) symmetry (set R q , see Table |T|). Results are shown in 
Fig. [|. Again, at high densities the Groningen and DD parameter set are much softer 
than the Bonn A parameter set. But their properties at low and intermediate densities 
are considerably different. The Groningen EoS exhibits an extreme softening at low den- 
sities. This is explained by the early appearance of the £~. The strong coupling to the S 
field, being only present in the Groningen parameterization, reduces the S~ effective mass 
= Ms- — To-s-^o- — 2r < 5 S -$,5. In contrast the isoscalar A hyperon does not couple 
to the isovector 5 meson. Moreover, because of isospin, the effective masses of the S + and 
5° are increased and their treshholds are shifted to higher density regions. This behavior 
is different from || where the inclusion of a constant coupling 6 meson did not noticeably 
change the EoS. 

Results for the composition of neutron star matter for the different parameter sets are 
presented in Fig. One recognizes that the early and strong appearance of the E~ temporar- 
ily reduces the electron fraction and strongly suppresses the myons. Clearly, this behavior 
opposes other results of RMF and Brueckner calculations, e.g. |5|,[/]]. The extraction of the 
isovector couplings from the Groningen potential leads to relatively large values indicating 
a slight increase of the 5 coupling above saturation density. This might cause the observed 
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behavior and could be fixed by adjusting the scaling factor RgT,- to the strength of the S 
isovector potential. Since experimental data is not available we used Rj:$ = R p y, = 2 from 
SU(6) symmetry instead. Further investigations are necessary. 

The composition of neutron star matter for the DD parameter set resembles the outcome 
of non-relativistic Brueckner calculations |)],[7|] with the S~ appearing before the A. For the 
Bonn A parameter set the A and £~ appear at relatively high densities of about 2.5 — 3po 
which also explains the stiffer equation of state. At higher densities when the other hyperons 
are present the number of hyperons exceeds the number of neutrons. This behavior is 
observed for all parameter sets with the strongest effect for the Groningen parameter set. 
The same result was also found in other models, e.g. |2||J. 

As mentioned in Sec. | one problem of RMF calculations with coupling constants fitted 
to properties of finite nuclei is the appearance of negative effective masses in neutron star 
matter at densities of p ~ 0.5 — 1 fm -3 ||. This is mainly caused by the small effective 
mass at saturation density that is needed in order to reproduce the observed spin-orbit 
splitting in finite nuclei This causes a rather stiff equation of state and confirms that 
the extrapolation of phenomenological parameter sets to high densities can be problematic. 
In the DDRH approach negative effective masses appear only at much higher densities where 
the validity of the model is already questionable because baryons and mesons will cease to 
be the relevant degrees of freedom. We therefore stop the calculation when negative effective 
masses appear and do not attempt to extrapolate to higher densities as in |]5|j22|. For the 
Bonn A parameter set m* gets negative around p ~ 1 — 1.2 fm -3 m 6 — 8po being mainly 
caused by the mentioned increase of the a coupling strength at high densities. For the 
Groningen parameter set the neutron effective mass gets negative at around p « 1.2 fm~ 3 . 
The reason is here the additional reduction of m* n by the coupling to the 5 meson. For the 
DD parameter set m* is still positive for the highest considered densities of p = 1.5 fm~ 3 . 
The different behavior of the DDRH model compared to constant coupling models originates 
in the density dependence of the vertices and is mainly caused by the fact that the vertex 
density dependence takes into account information from wider density range as found in 
finite nuclei and allows a safer extrapolation to higher densities. The medium dependence 
of the vertices implicitly includes higher order correction terms that can be partly identified 
with the nonlinear meson self-interactions of RMF models. Also, this explains why models 



with nonlinear vector self- interactions pT9|j20[| follow more closely DB calculations avoiding 
negative effective masses. 

Finally, we turn to the discussion of the results obtained with model 2. In Fig. [7] the 
equation of state for the different density dependent interactions is shown and Fig. || displays 
the corresponding particle distributions. Compared to model 1 the hyperons appear at very 
high densities of about 3po- I n addition their ratio increases very fast and around 3 — 4po 
the number of As already exceeds the number of neutrons. This leads to a stiffer EoS at low 
densities and a strong softening above e = 500 MeVfm -3 . The suppression of the hyperons 
in model 2 is a direct consequence of the different choice of the density dependence. The 
different rearrangement self-energies for each baryon multiplet induce a relative shift in the 
chemical potentials that impedes the appearance of the hyperons. The differences between 
the choices of scaling factors R a \ are the least pronounced for the Bonn A parameter set 
having the strongest density dependence and the most pronounced for the DD parameter 
set having the weakest density dependence. Here, the A does not appear for the densities 
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considered if the scaling factor set R q is chosen. This leads to the stiffest equation of 
state. For the semi-microscopic set R m all interactions have very similar equations of state 
approaching asymptotically the results of model 1 at high densities. Using set R q from 
SU(6) symmetry and including all octet hyperons we found that the appearance of the S 
and the 5 is shifted to even higher densities (p > 0.7 fm -3 ) or is completely suppressed for 
all interactions and does not noticeably modify the EoS. 

We conclude that the composition of neutron star matter strongly depends on the density 
dependence of the interaction. In practice this is problematic in model 2 since the low 
density behavior of the vertices has to be extrapolated from DB calculations and introduces 
additional uncertainties. On the other hand, there are indications that model 2 resembles 
the outcome of Brueckner calculations. Improvements of DB calculations in the low density 
region might resolve these problems. Again, we remark that we expect that a realistic density 
dependence is a mixture of model 1 and 2 and has to be determined from DB calculations 
of strange matter. 



IV. NEUTRON STARS 

The condition of hydrodynamical equilibrium inside a neutron star determines the con- 
nection between experimentally observable properties like the mass and possibly the radius 
of a neutron star and the theoretical equation of state of neutron star matter. In this work 
we only consider static spherically symmetric neutron stars. The condition of hydrodynami- 
cal equilibrium in general relativity is expressed by the Tolman-Oppenheimer-Volkoff (TOV) 
differential equations 



^P- = ^r\(r) (37) 
dr 

dP(r) [e(r) + P(r)] [m(r) + 47rP(r)r 3 ] 

- , (38) 



dr r 2 



i 2m(r) 



where r is the radial distance from the origin of the neutron star and m{r) is the mass 
contained in a sphere of radius r inside the star. The gravitational constant G has been set to 
1 for simplicity. Pressure P(r), energy e(r) and mass m(r) distributions can be calculated for 



a given pressure-energy relation P(e) as discussed in section 111 . The differential equations 
are solved by integrating from the center, starting with a central pressure P c , until the 
surface (radius r = R, P(R) = 0) of the neutron star is reached. The radius R of the star 
and its gravitational mass are related by M = m(R) = 4ir drr 2 e(r). Varying the central 
pressure one finds a mass-radius relation M(R) describing a family of neutron stars that 
will depend on the choice of the equation of state. 

Since it is a well known fact that Brueckner theory breaks down at very low densities 
one cannot rely solely on an EoS derived from the NN interaction for the description of 
neutron stars. Furthermore, in the surface region of neutron stars a so called crust of sub- 
nuclear density (where p is much lower than the normal nuclear matter saturation density 
Po) exists consisting of atoms and nuclei and having its own specific degrees of freedom. 
Even though the crust contributes only about 1% to the total mass and does not strongly 
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affect the maximum mass of the star, its thickness is about 10% of the radius and influences 
the size of very light stars. 

For the above reasons, we use the Baym-Pethick-Sutherland (BPS) EoS p3] for sub- 
nuclear densities p < 0.001 fm~ 3 and results of Negele and Vautherin for densities of p > 
0.001 fm -3 |p3| . The transition density to the DDRH EoS is defined by the intersection of 
both EoS. One finds values between 0.3 po and 0.7po depending on the different interactions. 

In Fig. [] mass-radius relations calculated with the equation of state of the Groningen 
parameter set are shown. For model 1 we find maximum masses of M max = 1.69M for the 
semi-microscopic scaling factor set R m and M max = 1.89Af Q for the SU(6) set R q . Including 
the full set of octet hyperons into the EoS the mass is further reduced to M max = 1.65M . 
The maximum density in the center of the star is similar in all three models and reaches 
p c ~ 0.91 — 0.99 fm~ 3 which corresponds to 5 — 6po- Including all hyperons the minimum 
radius of the neutron star decreases from -R m i n ~ 11.7 km to -R m i n — 10.8 km and the radius 
of neutron stars with the typical experimentally observed mass of M = 1.4M is reduced 
from R = 12.7 km to R = 11.7 km. The smaller radii are caused by the early occurrence of 
the E~ and the softening of the EoS at low densities that was discussed in Sec. |111 C. Also, 
the star mainly consists of hyperons as can be seen from the particle composition of Fig. |6|. 
Model 2 leads to significantly higher maximum masses of M max = 1.95M (R m ) respectively 
A^max = 2.12M Q (R q ) because of the appearance of the hyperons only at higher densities. 
Since, in addition, the maximum central density is only about p c ~ 0.73 fm -3 the star is 
mainly composed of nucleons. 

From these results the density dependence of model 1 seems to be favored because of the 
closer agreement with other calculations of neutron star masses. The big difference between 
the maximum masses for the two sets R q resp. R m emphasizes the sensitivity of the model 
to the choice of the scaling factors R a y- The effect is as large as the inclusion of addi- 
tional hyperons and should also be examined in RMF calculations with density independent 
coupling constants. 

For the Bonn A parameter set we find comparable results as is seen from Fig. [L0[ The 
maximum masses for model 1 of M max = 1.76M (R m ) resp. M max = 2.07M o (R q ) (only A) 
and the central density p c — 1.15 fm -3 are higher than for the Groningen parameter set due 
to the stiffer equation of state. Including the full set of octet hyperons the maximum mass 
is not reached before the effective mass gets negative at about p ~ 1.02 fm -3 where a mass 
of M = 1.81M is found. The same is observed for model 2 where the effective mass gets 
negative at about the same density indicating maximum masses higher than 1.7M (R m ) 
resp. 1.8M (R q ). Despite these differences the radii of neutron stars with 1.4M are all in 
the range of R = 12.5 — 12.8 km and in close agreement with the results for the Groningen 
parameter set. 



Results for the DD parameter set are displayed in Fig. 11]. Having the softest equation 



of state at high densities, this interaction also exhibits the smallest maximum masses of 
M max = 1.43M (R m ) resp. M max = 1.66M (R q ) (Model 1, only A) and M max = 1.44M 
(Model 1, all hyperons). For model 2 the masses are again higher (M max = 1.77M resp. 
M max = 2.O8M ) since the stars are mainly composed of nucleons. The DD parameter set 
reaches central densities up to p c ~ 1.3 fm -3 in model 1 and p c ~ 1.0 fm -3 in model 2. The 
soft equation of state reduces also the radius of the neutron stars. One finds a minimum 
radius of -R m i n — 10.3 km and a radius of R ~ 10.8 km for a mass of 1.4M if all hyperons 
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are included. 

Even though the three different parameter sets and the different choices of the density 
dependence give different predictions for the maximum mass of a neutron star the radius 
of a typical neutron star with a mass of M = 1.4M is in all models expected to be about 
12.5 km while the minimum radius can get as small as 10.3 km. This is explained by the 
fact that the radii are mainly influenced by the low to intermediate density behavior of the 
EoS which is similar in all considered models. On the other hand, the maximum masses are 
strongly influenced by the high density behavior being quite different for the investigated 
interactions and models 1 and 2. Our results are compatible with non-relativistic Brueckner 
calculations and relativistic Hartree-Fock calculations (fitted to the Bonn A NN 

potential) |22| that predict radii of 10 — 12 km and also possess a rather soft equation of state. 
In contrast, most relativistic mean- field models with density independent coupling constants 
find radii of 13 — 15 km, e.g. fl2",|55"|. Our maximum neutron star masses of 1.44 — 1.8M agree 



with the outcome of different RMF calculations while non-relativistic Brueckner calculations 
favor values of ~ 1.3M . This indicates that the non- relativistic EoS might be too soft 
at higher densities where relativistic effects become more important. The high value of 
~ 1.8M found for the Bonn A potential is probably caused by the inaccuracies from 
the polynomial extrapolation of the density dependent vertices to the high density region. 
Therefore, we favor the results of the Groningen and DD parameter set. We also performed 
calculations with a purely hadronic EoS to examine the total effect of the hyperons on the 
EoS. For the DD interaction a maximum mass of 2.O5M was found while the Groningen 
and Bonn A parameter sets have masses of 2.35M resp. 2.45M . The observed reduction 
AM ~ 0.5 — O.7M of the maximum mass caused by the inclusion of the hyperons is 
comparable to other models ]l]||[7j • 



V. SUMMARY AND CONCLUSION 



We have extended the DDRH model to the strangeness sector including the full set 
of SU(3)/ octet baryons. Using realistic density dependent interactions derived from DB 
calculations we examined the properties of neutron star matter. Starting with a medium- 
dependent NN interaction that is parameterized by density dependent meson-nucleon ver- 
tices we found that the extension of DB results from nuclear matter to hypermatter is 
not unique. The structure of the DB interaction strongly indicates that the ratio R a y of 
the nucleon and hyperon in- medium vertices and self-energies is mainly determined by 
the ratio of the corresponding free-space coupling constants and is only weakly affected by 
the background medium. We developed two different models to describe the connection 
of the density dependence of the baryon vertices to the surrounding medium. In model 1 
the vertices depend on the total baryon density and are not influenced by the strangeness 
asymmetry of the medium. This ensures that the relative strength of the vertices remains 
constant and follows closely the behavior of RMF calculations with constant couplings. An 
extreme assumption in this model is that the vertices are influenced equally from all baryons. 
However, in first order hyperons and nucleons should be independent of each other. This 
was investigated in model 2 that assumes a dependence of the vertices only on the baryons 
within the corresponding SU(3)/ multiplet. While this leads to a more realistic description 
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of slightly asymmetric systems, the model becomes problematic at very high asymmetries 
since it completely neglects the higher order effects from of the other baryon multiplets. In 
addition, in such systems the low density behavior turns out to be important which is not 
well determined from Brueckner calculations. We examined the properties of NA matter and 
found differences between the two models mainly at low densities and strangeness fractions 
fs that are nearly vanishing at high densities and values of fs close to one. The conclusion 
is that a realistic medium dependence of the in- medium vertices should be a mixture of 
both models including the asymmetry and the total density. Such a description is pending 
until DB calculations for the full strangeness octet will be available. 

Application of the DDRH theory to weak /5-stable matter leads to more complex con- 
dition for chemical equilibrium than in mean-field calculations with constant coupling con- 
stants. The occurrence of rearrangement contributions to the chemical potential connects 
all particle densities strongly to each other which complicates the practical calculation con- 
siderably. 

An important aspect in our calculations is the value of the scaling factors R a Y- Micro- 
scopic calculations of the cxA vertex show a strong deviation from the quark model value. 
Choosing either _R ct a = 0.49 from NA T-matrix results (scaling factor set R m ) or R^k = 2/3 
from SU(6) symmetry considerations (set R q ) and fixing the remaining scaling factors to the 
same experimental potential depth of the A we found strong differences in the equation of 
state. The smaller microscopic values lead to a much softer EoS and reduce the calculated 
maximum neutron star masses by 10-15%. This behavior even though enhanced by the 
density dependence of the interaction will also be visible for parameter sets with density 
independent coupling constants and should be closely examined. Microscopic calculations 
for the S and H vertices are not available restricting full scale neutron star calculations to 
the SU(6) symmetry values. But one has to review if these values are realistic and study 
effects of their modification on neutron star properties [|56| , |5"7ll . 

We remark that our calculations did not examine the influence of the hidden- strangeness 
mesons a s and 4> explicitly that are assumed to cause a highly attractive hyperon-hyperon 
interaction at intermediate densities and might lead to a small reduction of the neutron 
star radii. On the other hand, the couplings of these mesons posses some uncertainties that 
cannot be reliably fixed since experimental data are missing. In addition they are in part 
implicitly included in the microscopic scaling factors. 

An examination of neutron star properties favors model 1 as choice for highly asymmetric 
dense systems. While model 2 predicts too high maximum masses and a suppression of the 
hyperons, the results for model 1 are in close agreement with other calculations. For all 
examined density dependent interactions we find radii of 10 — 13 km for neutron stars 
with masses above 1.4M . This is also observed in non-relativistic Brueckner calculations 
while RMF calculations with phenomenological interactions usually favor larger radii. Our 
maximum masses of M max = 1.65M© for the Groningen parameter set and M max = 1.44M Q 
for the DD parameter set are confirmed by other RMF calculations and are slightly higher 
than in Brueckner calculations. 

We conclude that the DDRH model allows a consistent calculation of strange matter 
and neutron stars and yields results that are comparable with other models. It incorporates 
the properties of the DB model using microscopic interactions at various densities as input. 
The extrapolation to higher densities is more constrained than for phenomenological RMF 
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calculations that use only information from the limited density range of finite nuclei for the 
determination of their parameters. Improvements of the results, a more realistic density 
dependence of the hyperon-meson vertices and more restricted predictions will be possible 
if results from upcoming DB calculations for the SU(3) / baryon octet are available. 
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FIG. 1. Equation of state of symmetric nuclear matter for A admixtures from 0% to 40%. 
Shown are results for the Bonn A NN potential for different choices of the density dependence 
(model 1 and 2) and the scaling factors (set R q and R m ) . 
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FIG. 2. Binding energy (top) and saturation density (bottom) as a function of the strangeness 
fraction fs for the different choices of the A vertex. Results are shown for the Bonn A (left) and 
the Groningen (right) parameter set. 
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FIG. 3. Binding energy of pure A matter as a function of the density for the different choices 
of the A vertex. Results are shown for the Bonn A (left) and the Groningen (right) parameter set. 
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FIG. 4. Equation of state for neutron star matter in the DDRH model for different density 
dependent interactions. Results are shown for model 1 including only the A and the different 
scaling factors. The upper line of each interaction corresponds to calculations with set R q and the 
lower one to calculations with set R m . For details see text. 
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energy density e [MeV fm" 3 ] 

FIG. 5. Equation of state for neutron star matter for model 1 including all hyperons. Results are 
shown for the phenomenological density dependence of Ref . jljj (solid line) , Groningen parameter 
set (dotted line) and Bonn A parameter set (dashed line). The SU(6) scaling factors from set R q 
are used. 
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FIG. 6. Composition of /J-stable matter in the DDRH model for the phenomenological den- 
sity dependence of Ref. [45| (top), Groningen parameter set (middle) and Bonn A parameter set 
(bottom). Results are shown for model 1 including all hyperons and scaling factor set R q . 
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FIG. 7. Same as Fig. f|, but results are shown for model 2. 
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FIG. 8. Same as Fig. |6|, but the density dependence is implemented with model 2. The scaling 
factor set R m was used. 
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FIG. 9. Neutron star mass as a function of radius in the DDRH model for the Groningen 
parameter set. Results for model 1 including only A (solid lines), model 1 including all hyperons 
(dotted lines) and model 2 (dashed lines) are compared. The upper line of each model corresponds 
to calculations with set R q the lower one to calculations with set R m . For details see text. 
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FIG. 10. Same as Fig. |9] but for the Bonn A parameter set. 
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TABLES 



Rq 



model 1 


Rak 


RloA 


RcK,H 


Ru)k,Y, 


R<tE 


RujB 


Bonn A 


0.49 


0.4935 


0.6225 


2/3 


0.3618 


1/3 


Groningen 


0.49 


0.5217 


0.6061 


2/3 


0.3343 


1/3 


DD phenom. 


0.49 


0.5100 


0.6170 


2/3 


0.3421 


1/3 




Rm 






Rq 






model 2 


RaA 


Ru)K 


Ra\,T, 


Ru)k,Y, 


RaS 


RujE 


Bonn A 


0.49 


0.5690 


0.5690 


2/3 


0.3040 


1/3 


Groningen 


0.49 


0.5387 


0.5921 


2/3 


0.3222 


1/3 


DD phenom. 


0.49 


0.5251 


0.6061 


2/3 


0.3287 


1/3 



TABLE I. Scaling factor sets R q and R m for the a and lo hyperon-meson vertices of model 1 
and 2. For details see text. 
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